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Cost-Optimal Switching Protection Strategy in Adaptive Networks 
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Abstract —In this paper, we study a model of network adap¬ 
tation mechanism to control spreading processes over switch¬ 
ing contact networks, called adaptive susceptible-infected- 
susceptible model. The edges in the network model are ran¬ 
domly removed or added depending on the risk of spread 
through them. By analyzing the joint evolution of the spreading 
dynamics “in the network” and the structural dynamics “of 
the network”, we derive conditions on the adaptation law to 
control the dynamics of the spread in the resulting switching 
network. In contrast with the results in the literature, we allow 
the initial topology of the network to be an arbitrary graph. 
Furthermore, assuming there is a cost associated to switching 
edges in the network, we propose an optimization framework 
to find the cost-optimal network adaptation law, i.e., the cost- 
optimal edge switching probabilities. Under certain conditions 
on the switching costs, we show that the optimal adaptation 
law can be found using convex optimization. We illustrate our 
results with numerical simulations. 

I. Introduction 

Accurate prediction and effective control of spreading 
dynamics over networks are relevant problems in epidemi¬ 
ology and public health, computer malware, or security 
of cyberphysical networks. Although we find many recent 
advances in the field of network epidemiology [1], there 
are still many open questions to transfer this knowledge to 
realistic epidemiological situations. One fundamental result 
in the mathematical analysis of spreading in networks is the 
close connection between the eigenstructure of the contact 
network and epidemic thresholds [2]-[4]. This result enabled 
the authors in [5]-[8] to propose a convex optimization 
framework to design the optimal distribution of pharmaceu¬ 
tical resources to control disease spread. This framework is 
specially adapted to static network structures in which the 
pattern of interconnections does not change over time. As 
we argue below, this may not be the case in many practical 
situations. 

Social distancing is one of the most important nonphar- 
maceutical approaches to control disease spread over human 
contact networks [9], [10]. Examples of social distancing 
are, for instance, isolation of patients, school closures, and 
avoidance of crowds. In spite of the obvious effect that such 
behavior have on the dynamics of the spread, there is a lack 
of studies about the role of social distancing in the spread of 
diseases over human contact networks. One of the reasons is 
that social distancing induces an adaptation of the network 
structure that depends on the state of the infection. Although 
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there are results in the literature about disease spreading over 
time-varying networks (see, e.g., [11]-[13]), these works are 
based on the assumption that the evolution of the network is 
independent of the state of the individuals. In this paper, we 
propose a tractable framework to analyze the co-evolution of 
the state-dependent network structure and the dynamics of 
the spreading process taking place on it. 

Most of the available studies of spreading processes over 
human networks with social distancing have been relying 
on various unrealistic simplifying assumptions. The authors 
in [14]-[17] propose epidemic thresholds under the so- 
called mixing assumption; all the individuals in a network 
interact randomly with each other. However, this assumption 
is not satisfied in structured human populations. Although 
the analysis in [18] does not rely on the mixing assumption, 
it relies on the quantity called a reproduction number, whose 
validity for disease spread over time-varying networks is not 
yet fully established [19]. 

This paper analyzes, without the mixing assumption, the 
dynamics of spreading processes taking place in switching 
networks whose structure adapt to the state of the spread. The 
disease spread is modeled by an extended version of the well- 
known susceptible-infected-susceptible (SIS) model, which 
is called the adaptive SIS model [16]. Without the mixing 
assumption employed in [16], we derive conditions under 
which the network adaptation is able to protect against the 
spread of the disease. We furthermore use these conditions 
to propose a cost-optimal adaptation policy to contain the 
disease. This policy is based on the assumption that adapting 
the network structure to the state of the disease has an 
associated cost. The optimal policy can be then found by 
solving an optimization program. Under certain conditions, 
this optimization program can be effectively solved using 
elements from convex optimization [20]. 

This paper is organized as follows. In Section |I^ we 
introduce the adaptive SIS model studied in this paper. 
In Section III we analyze the exponential stability of the 
infection-free equilibrium of the adaptive SIS models. Based 
on our stability analysis. Sections IV and [V] study an cost- 
optimal adaptation strategy for networks of homogeneous 
and heterogeneous agents, respectively. 


A. Mathematical Preliminaries 

The probability of an event is denoted by P( ). The 
expectation of a random variable is denoted by £■[•]. We 
let I denote the identity matrix and Ip the p-dimensional 
vector whose entries are all one (we omit the dimension p 
when it is obvious from the context). A real matrix A, 
or a vector as its special case, is said to be nonnegative. 




denoted by A > 0, if all the entries of A are nonnegative. 
The notations A > 0, A < 0 and A < 0 are understood in 
the obvious manner. For another matrix B having the same 
dimensions as A, the notation A <B implies A—B<Q. We 
again understand A < B, A > B, and A > B in the obvious 
manner. The Kronecker product [21] of A and B is denoted 
by A (g) B. Let A be a square matrix. The maximum real part 
of the eigenvalues of A is denoted by We say that A is 

Hurwitz stable if rji^) < 0- Also, we say that A is Metzler if 
the off-diagonal entries of A are all non-negative. We say 
that A is irreducible if no similarity transformation by a 
permutation matrix makes A into a block upper triangular 
matrix. For matrices Ai, ..., A„, the direct sum 0;LiA,' 
is defined as the block diagonal matrix having the block 
diagonals Ai, ..., A„. When Ai, ..., A„ have the same 
number of columns, we define coli<,<„A, = col(Ai,... ,A„) 
as the block matrix obtained by stacking the matrices Ai, 
..., A„. 

A directed graph is a pair Q = (V,f), where V is a finite 
set of nodes, and f C V x V is a set of directed edges. Unless 
otherwise stated, we assume V = {1,.A directed path 
from i to j in Q is an ordered set of nodes (Iq, • • • ,ii) such 
that io = i, (4)4+i) G £ for k = 0, 1, and 4 = j- We 

say that Q is strongly connected if there exists a directed 
path from i to j for all i,j G V. The adjacency matrix of Q 
is defined as the n xn matrix A = such that aij = 1 

if {i,j) G £ and = 0 otherwise. Similarly, an undirected 
graph is a pair Q = (V,f), where V is a finite set and £ is 
a subset of unordered pairs {i,j} of the elements i,j G V. 
The adjacency matrix of an undirected graph is defined in a 
similar manner. A graph is strongly connected if and only if 
its adjacency matrix is irreducible. 

Finally, we recall basic facts about a class of optimization 
problems called geometric programs [20]. Let xi, ..., 
denote m real positive variables. We say that a real-valued 
function / of x = (xi,... ,Xm) is a monomial function if there 
exist c > 0 and ai,... G K such that /(x) = cxj' • • -x^"'. 
Also, we say that / is a posynomial function if it is a sum 
of monomial functions of x. Given posynomial functions /o, 
■ ■■, fp and monomial functions gi, g^, the optimization 
problem 

minimize /o(x) 

subject to fi{x) <1, / = 1,... ,p, 

8j{x) = l, j=l,...,q, 

is called a geometric program. It is known [20] that a geo¬ 
metric program can be converted into a convex optimization 
problem. 

11. Susceptible-Infected-Susceptible Model 
OVER Adaptive Networks 

This section introduces the model of spreading processes 
over adaptive networks studied in this paper and state the op¬ 
timal design problem under consideration. Each node in the 
network can be in one of two states: susceptible or infected. 
The state of node i evolves over time and is represented by a 


binary variable x; G {0,1}. We say that node i is susceptible 
at time t if x,(f) = 0, and is infected at time t if x,(t) = 1. In 
this paper, we model the evolution of x, as a continuous-time 
stochastic process taking values in {0,1}. We also assume 
that the structure of the network in which the spreading 
process is taking place evolves over time. In particular, we 
model the network ^ as a continuous-time random graph 
process taking values in the set of undirected graphs with n 
nodes. In other words, we model the dynamics of spreading 
as a stochastic process taking place over a random graph 
process. We denote by A/;(f) the set of neighbors of node i 
in the graph Q{t), i.e., A/}(f) = {j gV: {i,j} G Gf)}, and 
by A(f) = [aij{t)]ij the adjacency matrix of G{t). 

The spreading models over adaptive networks studied in 
this paper are formally introduced as the class of pairs 
{x,G) = i{xi}'l^i,G) satisfying the following definition: 

Definition 2.1: Let Go = (V,fo) be an undirected graph 
with adjacency matrix Ao = [aij{Q)]ij. The pair {x,G) is said 
to be an adaptive susceptible-infected-susceptible model over 
Go {ASIS model for short) if there exist nonnegative numbers 
Pi, 5i, (j)i, and \j/ij {i,j = such that the following 

conditions hold: 

a) G(0) = Go', 

b) The process (x,G) is Markov; 

c) For every i, the transition probabilities of x, are given by 

P{xi{t + h) = l\xi{t)=0)=Pi Y, Xkit)h-\-o{h), (1) 

kcMif) 

P{xi{t -\-h) = 0 I Xi{t) = 1) = 5ih-\-o{h), (2) 

where o{h) is a function such that \imk^Q o{h)/h = 0. 

d) For all /, j, the transition probabilities of aij are given by 

P{aij{t+h) = 0\aij{t) = l) = {<^iXi{t) + <^jXj{t))h + o{h), (3) 
P{aij{t +h) = \\aij{t)=0) = aij{0)^rijhPo{h). (4) 

e) yfij = xj/ji for all i and j. 

The constants jS,-, 5,, (j)i, and y/jj are respectively called 
infection, recovery, cutting, and rewiring rates. 

We can interpret the above model as follows. Item b) 
indicates that the future evolution of the spread, given the 
present state, does not depend on the past. The probabilities 
in c) describe how nodal states evolve. Notice that, if Git) 
were a static network, these probabilities would coincide 
with those of the NIMFA model [4] with heterogeneous 
infection and recovery rates. Eqn. ([T]) indicates that, if node i 
is susceptible and its neighbor j is infected, then i becomes 
infected with the instantaneous infection rate pi. Moreover, 
the rate is proportional to the number of infected neighbors. 
Eqn. (|^ implies that, once node i becomes infected, it will 
become susceptible with an instantaneous recovery rate 5,-. 

Item d) describes an adaptation mechanism of the net¬ 
work to the state of the disease. Eqn. 01 indicates that, 
whenever a node i is infected, the node adaptively removes 
edges connecting the node and its neighbors according to a 
Poisson process with rate 0,. This mechanism is designed 
to contain the spread through edges connected to infected 
nodes. Moreover, Q describes a mechanism for which ‘cut’ 




Fig. 1; Adaptively switching network. 


edges are ‘rewired’ or added back to the network. We assume 
that edge {i,j} is added to the network with a rewiring 
rate y/jj. See Fig. [T] for a schematic picture of these transition 
probabilities. Finally, Item e) follows from the assumption 
that Q{t) is undirected, although this is not an essential 
restriction and could be relaxed to account for directed 
contact networks. Also, notice that we have included the 
term a,;( 0 ) in 0 to guarantee that only those edges that 
were present at the initial time t = 0 can be added later on 
by the rewiring process. 

Remark 2.2: A model similar to the ASIS model proposed 
in this paper was studied in [16], where it was assumed that 
the initial graph Qo is the complete graph. A major difference 
between our model and the one in [16] is the information 
available to each node. In the model in [16], it is assumed 
that nodes know the states of their neighbors. In contrast, 
we do not assume to have access to this knowledge in our 
model. This difference has a direct implication in the link¬ 
breaking process. For example, in [16], an infected node does 
not break the edge between itself and its infected neighbors. 
On the other hand, in our model, an infected node will break 
edges independent of the state of its neighbors. 

Once the adaptive network under consideration is de¬ 
scribed, we define the exponential stability of the infection- 
free equilibrium pi{t) = 0 of ASIS models, as follows: 

Definition 2.3: For f > 0 let pi{t) = P{xi(t) = 1) be the 
infection probability of node i. We say that the infection- 
free equilibrium pfit) = 0 of the adaptive SIS model {x,Q) 
is exponentially stable if there exist K >0 and a > 0 such 
that p,(t) < Ke^“‘ for all /, f, and x,(0). We call a the decay 
rate. 

In many practical situations, there is a cost associated to 
the mechanisms of cutting and rewiring edges in a network. 
Accordingly, we assume we have two scalar cost functions / 
and g, defined on [ 0 ,°°), describing the cost associated to the 
rates of cutting and rewiring edges, respectively. The main 
purpose of this paper is to find the cost-optimal switching 
strategy, defined by the values of the cutting and rewiring 
rates, to drive the state of the spread towards the disease- 
free equilibrium at a given exponential rate. The total cost 
of a switching strategy is given by: 

C=£/(0,)+ Y. siVij)- 

1=1 {i,j}e£o 

We also assume the following bounds on the rates: 

W<¥ij<V (5) 


for some nonnegative numbers <j), (j), i//, and y/. Now, we are 


■0 


ready to state the problem investigated in this paper|^ 

Problem 2.4: Given a > 0, find the cutting and rewiring 
rates (j)i and ¥ij satisfying 0 such that the adaptive SIS 
model is exponentially stable with decay rate a and the total 
cost C is minimized. 

In this paper we solve Problem 2.4 under the following 
reasonable assumption: 

Assumption 2.5: Qq is strongly connected. Moreover, j3, > 
0 , 5i > 0 , and ¥ij > 0 for {Uj} G So. 

III. Stability Analysis 


In this section, we perform a stability analysis of the ASIS 
model {x,Q) over Qq. We begin by representing the model 
as a set of stochastic differential equations with Poisson 
counters. For 7 > 0, we let Ny denote a Poisson counter 
with rate 7 . We assume that all Poisson counters appearing 
in this paper are stochastically independent. We will use 
superscripts for the Poisson counters to distinguish those that 
has the same rates but are independent. Then, from Q and 
0 . the evolution of the nodal states can be described as: 

dxi = -Xi dNg. + ( 1 - X,) Y . ( 6 ) 

k€Ui(0) 

Similarly, from ([^ and 0 , the evolution of the edges can 
be written as: 

datj = -aij (x,- dN^^^ xj dN^^ ) + (1 - a/;) , (7) 

for all i and j such that {i,j} G So- 

Using the stochastic differential equations (|^ and Q, we 
derive an upper bounding linear model for the infection prob¬ 
abilities Pi. To state the linear model, we define the following 
variables. Let us define p{t) G ffi" by p = coli<K„p,'. Also, 
for i— and j € A//(0), define qij{t) = £’[a,j(f)x,(f)] 

and let qt — coljg_A7.(o) qtj and q = coli<K„q',. Let di denote 
the degree of node i in the initial graph Qq and m the 
number of the edges in Qq. Then, q has the dimension 
Yli=idi = 2m. We also introduce the following matrices. 
Define 7] G as the unique matrix satisfying: 

'Siq^ Y 4ki- (8) 

keJGiiO) 


Then define the matrices Bi = coli<,<„(]3,7]), 
B 2 = coli<K„(Ald, < 8 ) 7]), Di = ©"^1 5i, D 2 = ©Li(5,4)> 

^ = ©f=i(A4)> 'i’l = ©LiCol;eA/;.(o)r,> 

^2 = ©Li ®je5fi{0) ¥ij- Now, we can state the following 
theorem: 

Theorem 3.1: Define M G u^n+2m)xin+2m) 

M= ] (9) 

'Ll B2-D2-<P-'P2 


Then, for all xi(0), ..x„(0), it holds that 


p 

<M 

P 

.4. 


.4. 


( 10 ) 


* Since the design of infection and recovery rates have been previously 
studied in [5]-[8], we focus our attention on the design of (pi and only 
(although our framework can be easily extended to include /3, and 5, as 
additional design variables). 
















Proof: Taking the expectations in (|^ yields that 

[x/] = - 8iE[xi] +j5i Y, E[{1- Xi)aii,Xk] ■ 
keMi(0) 

Since ^[(l — Xi)aikXk] < E[aikXk] = qu, from the definition 
of Ti in ([U, we obtain dpi/dt < —5iPi + piTiq. This implies 
that dp/dt < —DipPBiq, which proves the upper half block 
of the inequality ( [T0| l. 

Then, let us evaluate dq/dt. The Ito rule for jump pro¬ 
cesses [22] yields that 

d{aijXi) = -aijXidN^^^ - aijXjXjdN^^^ + (1 - aifxidN^,^. 

- atjXi dNsi + aij (1 - x,) Y ‘^‘k^k . 

keMi(0) 

Taking expectations in this equation, we obtain 

^ = -^iE[aijXi] - (j)jE[aijXiXj] + 

at 

Xj/i )£■ [(1 - atj )x,] - 5iqij + Pi Y^ J (1 “ Xi)aikXk\. 

keMiiO) 


Since I.keAfiiO)E[aij{l - Xi)aikXk] < LkeAfiiO)E[aikXk] = fq, 
we obtain dqij/dt < ij/ijPi — {(j), + V^,;- -t- 5i)qij + PitjTiq from 
Stacking the variables qtj for all j € A/;(0) yields 
dqi/dt < colyg^.(o) {xj/ijPi) - ( 0 , + 5i)qi - Wjqi + PiiM ® Ei)q, 
where yfj = 0/gj V.(o) yfji- This proves the lower half block 
of the inequality ( | 10 | l and completes the proof. ■ 

From Theorem [m we immediately have the following 
sufficient condition for exponential stability of the infection- 
free equilibrium. 

Theorem 3.2: If M is Hurwitz stable, then the infection- 
free equilibrium of the adaptive SIS model is exponentially 
stable with a decay rate —T}{M). 

Before closing this section, we prove the following propo¬ 
sition that plays an important role in the rest of the paper. 

Proposition 3.3: The matrix M is irreducible. 

Proof: Define 


where 


L = 


O T 
J S 


7 = 01d,, T= col Ti, S= col (Id. 0 7]). (12) 

^ l<Kn KKn 

1=1 


Since Pi and y/ij are positive by Assumption |2.5[ if Mij = 0, 
then Lij = 0 for all distinct i and j. From this we see that, 
to show the irreducibility of M, it is sufficient to show that 
L is irreducible. 

In order to show that L is irreducible, we shall show 
that the directed graph Ti, defined as the graph having 
adjacency matrix L, is strongly connected. We identify the 
nodes -f 2 m of fL using the variables pi, ..., p„, 

qij (;GA/'i(0)), ..., qnj (;■ G A4(0)). Then, the upper- 
right block T of the matrix L shows that the graph fL has 
directed edge {pi,qji) for all and j G A/i(0). 

Similarly, from the matrices J and S, we see that 77 has 
the edges {qij,pi) and {qij,qki) for all i = and 


j,k G A/j(0). Then, let us show that fL has a directed path 
from Pi to pj for all i,j G Since is strongly 

connected, it has a path {io,...,ie) such that io = i and 
ie = j. Therefore, from the above fact, we can see that 
fL contains the directed path 

In the same way, we can show that 77 also contains the 
directed path {pi,qji,qij,Pi) for every {i,j} G So- These two 
observations show that 77 is strongly connected and, hence, 
L is irreducible. ■ 


IV. Homogeneous Case 


Based on the stability analysis presented in the previous 
section, we study the optimal design problem stated in 
Problem [2.4| We start our analysis by assuming that the ASIS 
model is homogeneous, as defined below (this restriction is 
relaxed in the next section); 

Definition 4.1: We say that the adaptive SIS model is 
homogeneous if there exist nonnegative constants j3, 5, (j), 
and yr such that pi = p, 5,- = 5, (pi = <j), and y/tj = yr for all 
i and j. 

In the homogeneous case, the stability criterion in Theo¬ 
rem reduces to the next simple condition. 

Theorem 4.2: Assume that the adaptive SIS model is 
homogeneous. Let p = rfi^o)- Then, the infection-free equi¬ 
librium of the adaptive SIS model is exponentially stable if 

PP-<I>-W I s/iPp + <l> + W)^-^P^ 


Proof: Assume that the model is homogeneous. Then, 
the matrix M defined in (|^ takes the form 


-51 PT 

y/J pS-{5 + <l) + y/)l\^ 


where the matrices J, T, and S are defined by ( [T2] i. We prove 
the theorem under the assumption that Pp f (p. Since Qo 
is strongly connected by Assumption |2.5| Aq is irreducible 
and therefore has a positive eigenvector v corresponding 
to the eigenvalue p (see [23]). Define the positive vec¬ 
tor w — coli<,<„(v/ld.). Then, the definition of 7] shows 
TiW = LkeMM = LkeMiiO) Vk = {Av)i = pvi and therefore 
Tw — Xv. In the same manner, we can show Sw = pw. Since 
we have Jv = w, for a nonnegative number c it follows that 


cv 


{Pp-c5)v 

w 


{cyr + pp - {5 + <p + yf))w 


(14) 


Hence, if a real number A satisfies the following equations: 


Pp -c5 = cX, cy/ +Pp - {5 + ^+ y) = (15) 

then, by ([T^, we see that the nonnegative vector col(cv,w) 
is an eigenvector of the irreducible and Metzler matrix M 
corresponding to the eigenvalue A. This implies that ri{M) = 
X (see [23, Theorem 17]). Therefore, the condition A < 0 is 
sufficient for exponential stability of the adaptive SIS model 
by Theorem 1 3.2 1 

To find such A, we solve ( [T5] l with respect to A and obtain 
X^ {25 (p yf — pp)X 5{5 (p yf) — pp{5 yt) = 0. 













This equation is satisfied by A = A+, where 

. ^ l3p-25-(t>-\i/ + ViPp + (t> + ¥f-4l3p(t> 

+ 2 

Then, the pair (c,A) = (j3p/(A++ 5),A+) satisfies ([T5|. 
We remark that A+ + 5 is positive because of the initial 
assumption )3p 7 ^ 0. Therefore, c > 0 and hence the above 
argument shows that 


t 7 (M) = A+. 


(16) 


Therefore, by Theorem |3.2| the infection-free equilibrium of 
the adaptive SIS model is exponentially stable if A+ < 0 , 
which is equivalent to ( [T3] l. ■ 

Remark 4.3: In the special case when the network does 
not adapt to the prevalence of infection, i.e., when 0 = 0 , 
Proposition |4.2| recovers the well-known stability condition 
5 > j3p(Ao) for the SIS models over static networks [2], [4]. 


The following theorem provides a solution to Problem 2.4 
in the homogeneous case: 

Theorem 4.4: Assume that the adaptive SIS model is ho¬ 
mogeneous. Let 0 and xf/ be the solutions of the optimization 
problem: 


minimize n/( 0 )-l-mg(t//) 

<!>•¥ 

subject to 0 > (jSt] — 5 + 1) (<//(5 — a) + 1), (17) 

0 < 0<05 V ^ V ^ V- 

Then, the pair (0, i/r) gives the solution of Problem |2.4| 
Proof: It is sufficient to show that ri(M) < — a if and 
only if the ([T7| holds, but this easily follows from ([T 6 |). ■ 


V. Heterogeneous Case 


In this section, we extend our analysis to non- 
homogeneous adaptive SIS models. We will show that 
Problem |2.4| can be effectively solved under the following 
assumption. 

Assumption 5.1: 

1) The values of xj/jj are given for every {i,j} G foi 

2) There exist constants r > 0 and s such that the func¬ 
tion F: [r—0,r — 0]—s + f(r — x) is a posyn- 
omial function. 

In order to state the main result of this section, we will 
need the next proposition. 

Proposition 5.2: Let 5 = max,- 5, and define 5; = 5 — 5,-. 
Similarly, let xj/ = max/ yV/y and define xfij = xj/— xj/ij. 
Let 01, ..., 0„ be real numbers. Define the matrices 
Di = ©Li5,', D2 = eUiSild,-), ^ and 

‘P 2 = ©"=1 ©;■£+] ( 0 ) Vij- Define the nonnegative matrix 

jQ^\Di + Wl + rI Bi 

'Ll B2+D2 + ^ + '^2 ■ 


Then, for a given a > 0, the following statements are 
equivalent: 

• There exist 0i,...,0„ S [0,0] such that ri(M) <—a. 

• There exist 0i,..., 0„ G [r — 0, r — 0] such that 77 (M) + 
a < V 4 ) + 5 + r. 


Moreover, between {0,}]Li and {0,}[Li, there is a one-to-one 
correspondence given by the equation 


01 = r- 0 ;. 


(18) 


Proof: Assume that there exist 0i,...,0„ G [0,0] sat¬ 
isfying v\{M) <—a. Define 0, by ( fTS] !. Then we see that 
M = M+ (5 + r+ v/)/. This implies Tj(M) + a<5 + r+xff. 
We also have 0, G [r — 0,r — 0]. The other direction can be 
shown in the same way and, hence, its proof is omitted. ■ 

to a 


Using Proposition 5.2 we can reduce Problem 2.4 


geometric program under Assumption 5.1 as stated in the 
following theorem: 

Theorem 5.3: Let 0i, ..., 0„ be the solutions to the 
following geometric program: 


minimize ^ + ( 0 ,) 

0.T ,=1 


subject to {M + ct/) V < (5 + r+v/) V, 

V > 0, 

r—^<^i<r—(j). 


(19a) 

(19b) 

(19c) 

(19d) 


Then, {0,}]Li, defined in ( [T 8 ] l solve Problem 2.4 

Proof: By Proposition 5.2 Problem 2.4 is equivalent 
to the optimization problem 


minimize ^f{r—^i) 
h i=i 

subject to T] (M + al) < 5 + r + 1 )/, 

7— 0 < 0 i<r— 0 , 


( 20 ) 


after the change of variables ( fTS] ). Minimizing the objective 
function in this problem is equivalent to minimizing the one 
in by the definition of F, whose constant term s can be 
ignored in the optimization. Then, since M+al is irreducible 
by Proposit ion |3.3| we can replace the constraint ( |20l l into 
( |19b| ) and ( |19c| l in the same way as in [ 8 ] using Perron- 
Frobenius lemma. Also, by a similar argument as in [ 8 ], we 
can show that ( [T^ is a geometric program. This is because 
F is a posynomial and each entry of the matrix M + al is 
a posynomial in the variables 0i, ..., 0„. The details are 
omitted. ■ 

Remark 5.4: When 17/,7 are also design variables, the 


above argument reduces Problem 2.4 to a signomial program, 
which are (in general) hard to solve [ 20 ]. 

VI. Numerical Results 

We illustrate our results with a numerical example. Let Qq 
be the graph of a social network of n = 247 nodes and m = 
940 edges. The adjacency matrix of the graph has spectral 
radius p = 13.53. We assume that all nodes have identical 
recovery rate 5 = 0.1 and infection rate j3 = 5/ (1.I p) = 

does 


4.2 


6.720 X 10^^. Since 5/j3 = (Ll)p > p. Theorem 
not guarantee the stability of the infection-free equilibrium 
when 0 = 0 , i.e., when the network does not adapt. 

Let us design the cost-optimal cutting rates so that the 
spread stabilizes towards the disease-free equilibrium in the 
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Fig. 2: Cost-optimal cutting rates for stabilization 


adaptive network. We assume 0 = 0, 0 = 4j3, and xj/ij = j3 and 
use the following cost function in our numerical simulations: 

We have chosen this function since it is increasing and 
presents diminishing returns. Also, we have normalized it, 
so that /(0) = 0 and /(0) = 1, and fixed r = 20. Let the 
desired exponential decay rate be a = 0.005 and solve the 
geometric program in Theorem |4.4| to obtain the optimal 
cutting rates 0,. Fig. shows a scatter plot for the optimal 
rates, 0,-, versus the degrees of the nodes for all i G V. The 
resulting switching policy suggests that, in general, nodes 
with a larger degree should have higher cutting rates (as 
could be naturally expected). However, the relationship be¬ 
tween the optimal cutting rates and the degrees is not trivial. 
Alternatively, we have also studied the relationship between 
cutting rates and other network centrality measures and K- 
scores (though we omit these figures for space limitations). 
Our simulations do not show any trivial dependency between 
cutting rates and any of the measures considered. 

VII. Conclusion 

In this paper, we have studied the dynamics of spreading 
processes taking place in networks that adapt their structure 
depending on the state of the dynamics. Our model is 
based on a collection of stochastic differential equations 
with Poisson jumps that model the joint evolution of the 
states of the process taking place in the network, as well 
as the evolution of the network structure. To illustrate our 
framework, we have focused our attention in a popular model 
of spreading dynamics, the SIS model, and study it dynamics 
over adaptive, switched networks. For this particular model, 
we have derived conditions for the dynamics of the spread 
to converge towards the disease-free equilibrium. Using this 
stability result, we have then formulated an optimization 
program to find the cost-optimal adaptive strategy to achieve 
stability. We have also showed that this optimization program 
can be efficiently solved using geometric programming. A 
numerical example was included to illustrate our results. An 
interesting future work is to fully investigate the difference 


of information structures in our model and the one in [16] 
addressed in Remark 
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